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We report on recent work concerning the phenomenology of initial-state parton- 
kr effects in direct-photon production and related processes in hadron collisions. 
After a brief summary of a fcr-smearing model, we present a study of recent results 
on fixed-target and collider direct-photon production, using complementary data 
on diphoton and pion production that provide empirical guidance on the required 
amount of kx broadening. This approach provides a consistent description of 
the observed deviations of next-to-leading order QCD calculations relative to the 
inclusive direct-photon and 7r° data. We also comment on the implications of these 
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q , results for the extraction of the gluon distribution of the nucleon 

Introduction 

> , 

00 \ Direct-photon production has long been viewed as an ideal process for measuring the 

■ gluon distribution in the proton [2]. The quark-gluon Compton scattering subprocess 

{gq—fyq) provides a large contribution to inclusive 7 production for which the cross 
sections have been calculated to next-to- leading order (NLO) [3] . The gluon distribution 
in the proton is relatively well-constrained at x < 0.1 by deep- inelastic scattering (DIS) 
• and Drell-Yan (DY) data, but less so at larger x. Consequently, direct-photon data 

' from fixed-target experiments can, in principle, provide an important constraint on the 

gluon content at moderate to large x. 

However, a pattern of deviations has been observed between the measured direct- 
photon cross sections and NLO calculations [4] . The discrepancy is particularly striking 
in the recently published higher-statistics data from E706 [5], for both direct-photon 
and 7T° cross sections. The final direct-photon results from UA6 [6] also exhibit evidence 
of similar, although smaller, discrepancies. The suspected origin of the disagreement is 
from effects of soft-gluon radiation. Such radiation generates transverse components of 



initial-state parton momenta, referred to in this discussion as fa?. (To be precise, hr 
denotes the magnitude of the effective transverse momentum vector, fcr, of each of the 
two colliding partons.) 

Evidence of significant kx has long been observed in production of muon, photon, 
and jet pairs. A collection of measurements of the average transverse momentum of the 
pairs {{pr)pair) is displayed in Fig. 1, for a wide range of center-of-mass energies (\/s). 
The values of (pT)pair are large, and increase slowly with increasing y/s. The values 
of (fcr) per parton (estimated as w (pr) pair / V2) indicated by these DY, diphoton, and 

1 This work was done in collaboration with L. Apanasevich, C. Balazs, M. Begel, C. Brombcrg, 
T. Ferbel, G. Ginther, J. Huston, S. Kuhlmann, A. Maul, J. Owens, P. Slattery and W. K. Tung [1]. 



© Institute of Experimental Physics, SAS, Kosice, Slovakia 



1 



2 



M. Ziclinski 



o 



D. 

V 




(GeV) 



Proton Data 



• Diphoton 
O Dimuon 
■ Dijet 



10, 



Vs (GeV) 

Fig. 1: (pr) of pairs of muons, photons, and jets produced in hadronic collisions versus ^fs. 



dijet data, as well as the inclusive direct-photon and 7r° production data, are too large 
to be interpreted as "intrinsic" — i.e., due only to the finite size of the proton. (From 
the observed data, one can infer that the average kr per parton is about 1 GeV/c 
at fixed-target energies, increasing to 3-4 GeV/c at the Tevatron collider, while one 
would expect {kr) values on the order of 0.3-0.5 GeV/c based solely on proton size.) 
Perturbative QCD (PQCD) corrections at NLO level are also insufficient to explain the 
size of the observed effects, and, in fact, full resummation calculations are required to 
describe DY and W/Z [7, 8, 9], and diphoton [10, ff] distributions. Similar soft-gluon 
(or kr) effects can be expected in all hard-scattering processes, such as the inclusive 
production of jets or direct photons [12, 13]. 

After reviewing a phenomenological model for kx effects in direct-photon production 
we will discuss its applications to data, as well as the implications for determining the 
gluon distribution. A more detailed presentation of these results can be found in [1]. 



kr Smearing Model 

The Collins-Soper-Sterman resummation formalism [14] provides a rigorous basis for 
theoretical understanding of soft-gluon radiation effects. Despite recent progress [15, 
16], no full treatment of inclusive direct-photon cross sections is yet available. In its 
absence, we use a PQCD-based model that incorporates transverse kinematics of initial- 
state partons to study the major consequences of kr for direct-photon production (and, 
by extension, for all hard-scattering processes). 
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In PQCD, the expression for the leading-order (LO) cross section for direct-photon 
production at large px has the form: 

<?(hih 2 -> jX) = j dx x dx 2 fa^htixi.Q 2 ) fa 2 /h 2 {x 2 ,Q 2 )^{a!a 2 -> 703), 

where a is the hard-scattering matrix element, and f ai /h! an d f a2 /h 2 are the parton 
distribution functions (pdf) for the colliding partons a\ and a 2 in hadrons hi and h 2 , 
respectively. To introduce kx degrees of freedom, we extend each integral over the 
parton distribution functions to the fcr-space, 

dx f a /h(x, Q 2 ) -> dxd 2 k T g(k T ) f a / h {x, Q 2 ), 
and take g(k T ) to be a Gaussian (as justified, e.g., by E706 data on high-mass pairs), 

9{kT) = ^(kjr- 

Here, (k T ) is the square of the 2-dimensional RMS width of the kx distribution for one 
parton and is related to the square of the average of the absolute value of kx of one 
parton through (k T ) = 4(k T ) 2 /7r. 

Since an exact treatment of the modified parton kinematics can be implemented in a 
Monte Carlo framework, but is more difficult in an analytic approach, it is convenient to 
employ Monte Carlo techniques in evaluating the cross sections according to the above 
prescription. In general, because the unmodified PQCD cross sections fall rapidly with 
increasing px, the net effect of kx smearing is to increase the expected yield. We denote 
the enhancement factor as K(px). 

A Monte Carlo program that includes such a treatment of kx smearing, and the 
LO cross section for high-p^ particle production, has long been available [17]. The 
program provides calculations of many experimental observables, and can be used for 
direct photons, jets, and for single high-py particles resulting from jet fragmentation 
(such as inclusive 7r° production). Unfortunately, no such program is available for NLO 
calculations, but one can approximate the effect of kr smearing by multiplying the NLO 
cross sections by the LO fcT-enhancement factor. Admittedly, this procedure involves 
a risk of double-counting since some of the /cr-enhancement may already be contained 
in the NLO calculation. However, we expect such double-counting effects to be small. 

A complete treatment of soft-gluon radiation in high-p^ production, should eventu- 
ally predict the effective kr values expected for each process and sfs. We will employ 
(kx) values representative of those found in comparisons of kinematic distributions in 
high-mass pair data with the above described model. 

Applications of the kr Model to Data 

The experimental consequences of kx smearing are expected to depend on the colli- 
sion energy. At the Tevatron collider, the smallest photon px values probed by the 
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CDF and DO experiments are rather large (10-15 GeV/c), and the fey-enhancement 
factors modify only the very lowest end of the px spectrum, where px is not signifi- 
cantly greater than fey. In the E706 energy range, large fey-effects can modify both the 
normalizations and the shapes of the cross sections as a function of px- Consequently, 
E706 data provide a particularly sensitive test of the fey model. At lower fixed-target 
energies, the fey enhancements are expected to have less px dependence over the range 
of available measurements, and can therefore be masked more easily by uncertainties in 
experimental normalizations and/or choices of theoretical scales. Nonetheless, the UA6 
and WA70 data generally support expectations from kx smearing. 

Comparisons to Tevatron Collider Data 

At the Tevatron collider, the above model of soft-gluon radiation leads to a relatively 
small modification of the NLO cross section. In Fig. 2 we compare the CDF and 
D0 isolated direct-photon cross sections [18] to theoretical NLO calculations with and 
without kx enhancement. In the lower part of the plot we display the quantity (Data- 
Theory)/Theory, for NLO theory without the fey-enhancement factor, and the expected 
effect from fey for (fey) = 3.5 GeV/c. This is the approximate value of (fey) per parton 
measured in diphoton production at the Tevatron [18], and one expects a similar (fey) 
per parton for single-photon production. (In the diphoton process, the 4- vectors of the 
photons can be measured precisely, providing a direct determination of the transverse 
momentum of the diphoton system, and thereby (fey).) 

As seen in Fig. 2, the fey effect diminishes rapidly with increasing px, and is essen- 
tially negligible above «30 GeV/c. The trend of deviations of NLO calculations from 
the measured inclusive cross sections is described reasonably well by the expected fey 
effect. Some of the observed excess can be attributed to the fragmentation effects in 
the isolated direct-photon production [19], but this alone cannot account for the entire 
deviation of the theory from data. 

Comparisons to E706 Data 

The conventional ((fey)=0) NLO calculations yield cross sections that are signficantly 
below the E706 direct-photon and 7T° measurements [5] (see Fig. 3). No choices of cur- 
rent parton distributions, or conventional PQCD scales, provide an adequate description 
of the data (for the presented comparisons all QCD scales have been set to py/2). The 
previously described fey-enhancement algorithm was used to incorporate the effects of 
soft-gluon radiation in the calculated yields. That is, the theory results plotted in the 
figures represent the NLO calculations multiplied by fey-enhancement factors K(px). 

As seen in Fig. 3, the NLO theory, when supplemented with appropriate fey en- 
hancements, is successful in describing both the shape and normalization of the E706 
direct-photon cross sections at both y/s — 31.6 GeV and 38.8 GeV. As expected, the 
fey-enhancement factors affect the normalization of the cross sections, as well as the 
shapes of the px distributions. The values of (fey) = 1.2 GeV/c at y/s = 31.6 GeV, 
and 1.3 GeV/c at y/s = 38.8 GeV, provide good representations of the incident-proton 
data. Both (fey) values are consistent with those emerging from a comparison of the 
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Fig. 2: Top: The CDF and D0 isolated direct-photon cross sections, compared to NLO theory 
without kr (dashed) and with kr enhancement for (kr) = 3.5 GeV/c (solid), as a function 
of pt- Bottom: The quantity (Data-Theory) /Theory (for theory without kr adjustment), 
overlaid with the expected effect from kr enhancement for (kr) = 3.5 GeV/c. 



same PQCD Monte Carlo with E706 data on the production of high-mass 7r°7r°, 771- , 
and 77 pairs [5] . Similar conclusions are reached from comparisons between calculations 
and E706 data for 7r~Be interactions at y/s = 31.1 GeV (not shown). 

For comparison, results of calculations using (kr) values ±0.2 GeV/c relative to the 
central values are also shown in Fig. 3. These can be taken as an indication of uncer- 
tainties on (kr), on the basis of several considerations. These include: (i) the range 
of kr values inferred from different distributions in the E706 high-mass pair data; (ii) 
differences observed between photon and tt° results; (iii) comparisons of (kr) values re- 
quired in inclusive cross sections with those representing the properties of massive pairs 
at E706 and WA70/UA6 energies, and (iv) the differences between dimuon, diphoton, 
and dijet values of (pr)pair seen in Fig. 1. Examples of enhancement factors K(px) for 
fixed-target experiments are displayed in Fig. 4. 

It is interesting to note that, for the fixed-target energy range, the kx enhancement 
increases at the highest values of pr- The shape of K{px) can be understood through 
the following argument. At the low pr end of the measurements, a (kr) of ml GeV/c 
is non-negligible in comparison to the px in the hard-scattering, and the addition of 
kx smearing therefore increases the size of the cross section (and steepens the slope). 
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Fig. 3: The photon and tt° cross sections from E706 at ^fs = 31.6 GeV (left) and 38.8 GeV 
(right) compared to /Vr-enhanced NLO calculations. 
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Fig. 4: Left: The variation of kr enhancements, K(pr), relevant to the E706 direct-photon 
data for protons at = 31.6 GeV, for different values of {kr}- K(pr) used by the MRST 
group [20] is also shown. Right: The same for y/s — 23.0 GeV (relevant to WA70 data). 



At highest pt (corresponding to large x), the unmodified NLO cross section becomes 
increasingly steep (due to the rapid fall in parton densities), and hence the effect of kx 
smearing again becomes larger. 

NLO calculations for 7r° production have a greater theoretical uncertainty than those 
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for direct-photon production since ir° production involves parton fragmentation. How- 
ever, the kx effects in 7r° production can be expected to be similar to those observed 
in direct-photon production, and the 7r° data can be used to extend tests of the conse- 
quences of kx smearing. Figure 3 also shows comparisons between NLO calculations [21] 
and 7T° production data from E706, using BKK fragmentation functions (ff) [22]. The 
previously described Monte Carlo program was employed to generate fc^-enhancement 
factors for 7r° cross sections, and (kr) per parton values similar to those that provided 
good agreement for direct-photon data also provide a reasonable description of n° data. 
For 7T° production, an additional smearing of the transverse momentum expected from 
jet fragmentation has also been taken into account. 



Comparisons to WA70 and UA6 Data 

Both WA70 and UA6 have measured direct-photon production with good statistics, 
and their data have been included in recent global fits to parton distributions. WA70 
measured direct-photon and 7r° production in pp and ir~p collisions at ^/s = 23.0 GeV 
[23], and UA6 has recently published [6] their final results (with substantially reduced 
uncertainties) for direct-photon production in pp and pp collisions at = 24.3 GeV. 
These center of mass energies are smaller than those of E706, and the (kr) values are 
therefore expected to be smaller (perhaps of order 0.7-0.9 GeV/c, based on Fig. 1). 
WA70 has compared kinematic distributions observed in diphoton events (for tt~p in- 
teractions) to NLO predictions, and has found that smearing the NLO theory with an 
additional (kx) of 0.9±0.2 GeV/c provides agreement with their data [24]. We therefore 
use this (kx) as the central value for the ^-enhancement factors for both experiments, 
and vary the (kx) by ±0.2 GeV/c, as was done with E706. Over the narrower px range 
of WA70 and UA6 measurements, the effect of kx is essentially to produce a shift in 
normalization, as illustrated on the right side of Fig. 4. 

Comparisons of the WA70 direct-photon and 7r° cross sections with the fc^-enhanced 
NLO calculations are shown in Fig. 5 (using QCD scales ofpx/2)- The tt° cross sections 
both for incident proton and n~ beams, and the photon data from incident n~ beam, 
all lie above the NLO calculations for (kx}— 0, and are in better agreement with the 
/cT-enhanced calculations; only the photon cross section for incident protons seems not 
to require a kx correction. 

The photon and tt° cross sections from UA6 for pp and pp scattering are shown in 
Fig. 6. The photon cross section for pp interactions lies clearly above the NLO calcula- 
tion for (kx) = 0, but is consistent with fcr-adjusted calculations for (kx) in the range 
of 0.7-0.9 GeV/c. The result for pp interactions is also above the unmodified NLO 
calculation, but requires a smaller value of (kx}- We note that the dominant produc- 
tion mechanisms for the two processes are different: quark- gluon Compton scattering 
dominates for pp, and qq annihilation for pp at the UA6 energy. As in the case of E706 
and WA70, the UA6 7r° cross sections are higher than the NLO calculation without kx, 
and can be described much better by introducing kx enhancement. 
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Fig. 5: The photon and 7r° cross sections from WA70 at yfs = 23.0 GeV for incident protons 
(left) and 7r~ (right), compared to fcT-enhanced NLO calculations. 




Fig. 6: The photon and 7r° cross sections from UA6 at y/s — 24.3 GeV for incident protons 
(left) and antiprotons (right), compared to fcT-enhanced NLO calculations. 



Impact on the Gluon Distribution 

It is now generally accepted that the uncertainty in the gluon distribution at large x is 
still quite large. Thus, it would appear important to incorporate further constraints on 
the gluon, especially from direct-photon data. To investigate the impact of k T effects 
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Fig. 7: A comparison of the CTEQ4M, MRST, and CTEQ4HJ gluons, and the gluon distri- 
bution derived from fits that use E706 data. The g\ and g[ gluon densities correspond to the 
maximum variation in (kr) that MRST allowed in their fits. 

on determinations of the gluon distribution, we have included the E706 direct-photon 
cross sections for incident protons, along with the DIS and DY data that were used in 
determining the CTEQ4M pdfs, in a global fit to the parton distribution functions. The 
CTEQ fitting program was employed to obtain these results [25], using the NLO PQCD 
calculations for direct-photon cross sections, adjusted by the fcr-enhancement factors. 
However, the WA70, UA6, CDF, and D0 data were excluded from this particular fit. 
The resulting gluon distribution, shown in Fig. 7, is similar to CTEQ4M, as might have 
been expected, since the fcr-enhanced NLO cross sections using CTEQ4M provide a 
reasonable description of the data shown in Fig. 3. 

The new MRST gluon distribution [20] (also shown in Fig. 7) is significantly lower 
than CTEQ4M at large x. While the MRST fit employs fcy enhancements (obtained us- 
ing an analytic integration technique) , it attempts to accommodate the WA70 incident- 
proton direct-photon data, which does not exhibit an obvious kr effect. In addition, 
the MRST /cr-enhancements are larger than ours at large pr (see Fig. 4), resulting 
in a smaller gluon at large x [26]. In contrast, the CTEQ4HJ gluon distribution [27], 
designed to improve the description of the high-py jet data from CDF in Run IA, is 
much larger than CTEQ4M in the same x range. This spread of the solutions for the 
gluon distribution at large x is uncomfortably large, and additional theoretical effort is 
warranted to properly incorporate the available direct-photon data in the pdf fits. 
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Conclusions 

We have described a phenomenological model in which (for) values used in the calcula- 
tions of for-enhancements are derived from data. The results are remarkably successful 
in reconciling the data and theoretical calculations for a broad range of energies. The 
for-enhancements improve the agreement of PQCD calculations with E706, UA6, and 
7T~ beam WA70 direct-photon cross sections over the full px range of measurements, 
as well as with the low-pr end of CDF and D0 results. All discussed fixed-target 7r° 
measurements also agree much better with such for-enhanced calculations. 

A definitive conclusion regarding the quantitative role of for effects in hard scattering 
awaits a more rigorous theoretical treatment of soft-gluon radiation. Such theoretical 
progress is crucial for a more reliable determination of the gluon distribution, especially 
in the large- a: region, where significant uncertainties remain. 
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